Photonic Band Gap in the Triangular Lattice of BEC vortices 
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We investigate the photonic bands of an atomic Bose-Einstein condensate with a triangular vortex 
lattice. Index contrast between the vortex cores and the bulk of the condensate is achieved through 
the enhancement of the index via atomic coherence. Frequency dependent dielectric function is used 
in the calculations of the bands, resulting in photonic band gap widths of a few MHz. 
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I. INTRODUCTION 

A rotating Bose-Einstein Condensate (BEC) manifests 
the formation of vortices after a critical rotation fre- 
quency. Furthermore, constituent vortices exhibit a pe- 
riodic structure, which is generally a triangular lattice 
[E HI H Q ■ Even near the borders of the condensate, lat- 
tice distortion is small [4(. From the theoretical point of 
view, a rapidly rotating BEC can be treated analytically 
and the density is found to be the product of a slowly 
varying function and a periodic function [H, 0] • 

Usual imaging of the vortices is carried on during the 
ballistical expansion of the condensate [l|, [H, H, H, 0] ■ An 
in-situ imaging was accomplished recently, by visualizing 
the 2D image of the lattice along the rotation axis, while 
condensate is in the trap Q . 

In the BEC experiments, rotation frequency is not di- 
rectly measured, but deduced from the change in the as- 
pect ratio of the cloud 0, H H E 0- 

However, we re- 
cently proposed a method, which is based on the reflec- 
tion through the directional pseudo photonic band gaps 
Q. Photonic band gaps in a triangular vortex lattice 
are obtained through the enhancement of the refractive 
index via quantum coherence [j^.[Tl| so that sufficient in- 
dex difference is generated between the vortex cores and 
the bulk of the condensate. Without an index enhance- 
ment scheme, a usual rotating BEC with a vortex lattice 
cannot exhibit high enough index contrast to obtain pho- 
tonic band gaps. BECs are rather dilute, and being an 
atomic gaseous medium, they exhibit dispersion only in 
highly absorptive regimes. When the ground and excited 
state are coupled to other auxiliary levels, however, ab- 
sorbtion in the resonant transition of the probe beam can 
be cancelled. This is due to the quantum interference of 
various absorbtion paths. This way, one can benefit from 
large dispersion at an atomic resonance, without absorp- 
tion. 

Utilizing the index enhancement scheme, triangular 
lattice of BEC vortices can generate both directional and 
complete photonic band gaps. Directional pseudo-gaps 
are also called as partial or stop gaps [l2| . Radiation can- 
not propagate in certain directions determined by these 
pseudo-gaps, but can in others. The rotation frequency 
of the condensate can be measured from the chopping in 



the reflected or transmitted probe beam at a directional 
photonic band gap. 

In a previous study, Q we have demonstrated the pres- 
ence of a photonic band gap within the frequency win- 
dow of index enhancement. In this paper, we discuss the 
photonic band structure for the full frequency regime, 
extending our work beyond the index enhancement win- 
dow. 

Though our general examination simply verifies exis- 
tence of the band gap within the index enhancement win- 
dow, beyond this region electric susceptibility becomes a 
complex valued function of frequency, for which definition 
of photonic band gaps are not straightforward. In the 
particular index enhancement scheme we consider here, 
there are absorption and gain regimes, where we have 
found no photonic band gaps. This allows for selective 
stoppage of the probe pulse among the other electromag- 
netic fields that are in use for trapping the condensate 
atoms and for index enhancement schemes. 

Despite the dense literature on photonic crystals, stud- 
ies of photonic energy bands which take complex, fre- 
quency dependent dielectric constants into account are 
rather sparse. Few recent studies investigate materials 
with small absorption [l3[ and low filling fractions of the 
dispersive and absorptive component 0, HEI. [lr| . More 
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recent works go beyond these limitations 
These efforts focus on understanding the properties of 
photonic crystals fabricated from metallic materials, 
which can have large complete photonic band gaps in the 
visible region of the electromagnetic spectrum. Drude- 
like models of metallic components within a dielectric 
host is used in the modelling of such crystals. Absorption 
is of negligible importance within transparency window 
of the Drude model which is about half of the plasma 
frequency. Small, yet realistic, amounts of absorption 
hardly changes the band structure . In the regions of 
appreciable absorption however, there may be no band 
gaps [l8j]. Our results, for a rotating atomic BEC, are 
in agreement with these results for photonic crystals of 
metallic materials. We note that, quite recently, broad- 
band absorptive properties of metallic photonic crystals 
are found to be advantageous for various applications 
[20| . Similarly, gain regime is important for understand- 
ing lasing properties of photonic crystals [2l| . Index en- 
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hancement schemes we suggest to use for rotating BEC, 
offer absorptive as weii as gain regimes beyond the index 
enhancement window, associated with lasing without 
inversion, 

The understanding of photonic crystals with complex 
dielectric constants is not trivial, due to the lack of well- 
defined group velocity for complex energy bands [22| . In 
this paper, we limit ourselves to the determination of the 
band structure, and do not discuss the details of beam 
propagation beyond the existence of the band gaps. 

The paper is arranged as follows. In section |TT] we 
overview the upper-level microwave scheme, which leads 
into the index enhancement with vanishing absorbtion. 
We introduce the system parameters and the resulting 
dielectric susceptibility. In section UTTT we obtain the ma- 
trix equations from the Master Equation of the photonic 
crystals and illustrate the method of solutions for the 
case of frequency dependent, complex dielectric function. 
In section. IIV1 we present the resulting photonic bands 
for two different lattice parameters and then discuss the 
properties of the photonic bands. Section[V]is a summary 
of our results. 



II. DIELECTRIC FUNCTION OF THE VORTEX 
LATTICE 

In this section, we describe the atomic coherence and 
path interference effects leading to high index of refrac- 
tion with vanishing absorbtion. We review the deriva- 
tions of Ref. [l(| in a compact form and describe the 
physics of the system. We calculate the real and imagi- 
nary parts of the dielectric susceptibility. 

There are various index enhancement schemes [Io| . 
Among them, we specifically consider upper-level mi- 
crowave scheme as it leads strong index contrast, though 
Raman scheme may also be useful due to its wider fre- 
quency window for index enhancement (Io| . The corre- 
sponding level diagram for upper-level microwave scheme 
is shown Fig. [TJ A weak optical probe field, E, of fre- 
quency u is coupled to the two levels a and b through 
electric dipole interaction. A third level c is also coupled 
to level a via a strong resonant microwave field of Rabi 
frequency f2 M . Coupling level a to level c allows for dif- 
ferent possible paths of absorption. Destructive quantum 
interference of these paths may cancel the absorption of 
the probe field at a certain frequency At the same 
time, a high refractive index can be generated by main- 
taining some population in level a to ensure high dipole 
moment for levels a and b. Indirect pump mechanisms 
are introduced to realize that. Parameters 77, r and r c 
are the pump rates from level a to c, b to a and b to c 
respectively. 

In Fig. [TJ 7^, 7 and j c denote the decay rates (inverse 
lifetime) of levels c to a, a to b and c to b respectively, due 
to the collisions and radiation. We can consider 7 M , 7 C <C 
7, as 7 C and 7 M are for dipole forbidden and microwave 
transitions, respectively. 




FIG. 1: Upper-level microwave scheme for index enhancement 
[Tol | . Upper two levels a and c are coupled via a strong mi- 
crowave field of Rabi frequency fi M . Weak probe field E, of 
optical frequency uj is coupled to levels a and b. Decay (7) 
and pump (r) rates are indicated. 



Following Ref. [T3|, we set 77 = r = 0, and choose r c = 
f2 M =7. We neglect 7^ and j c . In that case, frequency 
dependent electric susceptibility x{ u ) — x'( w ) + w"( w ) 
is a complex function of frequency, with its real x' an d 
imaginary parts x" are given by [l0( 
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where N is the number density of atoms and A is the 
wavelength of the optical transition a — > b. We define a 
dimensionless frequency 



(3) 



centered at the resonance frequency to a b and scaled with 
the decay rate of atomic coherence. 

Susceptibilities, given by Eqs. {TJ) and j2j), are for a 
dilute condensate. In the case of a dense condensate, 
first correction to the susceptibility is equivalent to a local 
field correction [23j in the form 



XlocO) 



i-xH/3 



(4) 



Real and imaginary parts of the corresponding dielectric 
function, 



= (nj) = 1 + XiooM 



(5) 



are plotted in Fig. [2] for a Rubidium 87 gas. The vertical 
line indicates the enhancement of the polarization at the 
frequency of vanishing absorption. Here, we define fio = 
2.37 x 10 15 Hz as the frequency at which absorption is 
zero. We also define wq = (f^o — ^ab)H — 1-22 as the 
corresponding value of the scaled and shifted frequency 
vo. We employ these definitions throughout the paper. 

Isotopes of alkali metals are typically used in the BEC 
experiments, and we specifically consider the energy lev- 
els of Rubidium. Fine-structure energy levels of Rubid- 
ium, which corresponds to b, a, and c levels of Fig. [TJ are 
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5si/2, 5pi/2, and 6si/ 2 , respectively. The wavelengths of 
a — b and a — c transitions become A = 794 nm and 
= 1.32 /im. The lifetime of the probe resonance 
level (Spin) is 27 ns which corresponds to the decay rate 
7 = 2?r x 6 MHz. 

In the vicinity of the center of the condensate cloud, 
the dielectric function plotted for the peak density in 
Fig. [5] can be assumed. When the vortices are present 
in this central region however, their spatial profile will 
influence the dielectric function. Assuming a dilute ther- 
mal gas background at ultracold temperatures, index en- 
hancement will be influential on the dense condensate 
only. Density of the condensate drops rapidly to zero 
at the vortex positions. Dielectric constant within the 
vortex core can be supposed to be same with that of 
vacuum, 6q — 1. Width of vortices is in the order of 
the coherence length of a condensate, which is given 
by £ = 1/ \/8TrNa sc , where a sc is the s-wave scattering 
length of the inter-atomic collisions. The typical val- 
ues of the coherence length is a few hundred nanometers 
and smaller than the optical wavelength. Spatial mod- 
ulations on the dielectric function can be introduced by 
( = (0 + (eioc — £ o)p(^) so that 

e l00 (f, w) = l + p(f)\ loc (w) . (6) 

Here p(r) stands for the normalized spatial profile of a 
vortex, and r is the radial distance from the center of the 
vortex core. This model dielectric function (|6|) drops to 
vacuum value at the cores of the vortices and recovers its 
bulk value in a few coherence lengths. 

In the band calculations for a triangular lattice of vor- 
tices, we considered a hexagonal Wigner-Seitz unit cell, 
which contains a single vortex core at the center of the 
unit cell. We used Pade's analytical form, derived in Ref. 
(24| , for the vortex density profile 

r 2 (0.3437 + 0.0286r 2 ) 

o(r) = — (7) 

h ' 1 + 0.3333r 2 + 0.0286r 4 ' v ' 

where r is scaled with the coherence length £. This den- 
sity behavior is valid in one unit cell. p(r) becomes zero 
at the center and goes to 1 towards the edges of the unit 
hexagonal cell. We choose the lattice constant a in terms 
of the coherence length to fix the filling factor of the vor- 
tices. We used two different values, a = 10£ and a = 4.5£, 
in the computations. 

In the following sections, we discuss the propagation 
of probe beam through the vortex lattice which has a 
dielectric function given by Eq. ([6]) and look for possi- 
ble band gap formations about the index enhancement 
frequency. 

III. CALCULATION OF THE PHOTONIC 
BANDS 

The stable lattice type for a single component rotating 
BEC is triangular [l|, 0] ■ Density profile is composed of 
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FIG. 2: Real (solid-line) and imaginary (doted-line) parts 
of local dielectric function ei oc (u;) as a function of scaled 
frequency w = (w — aj a b)/j, for the particle densities (a) 
N = 5.5 x 10 20 m -3 and (b) N = 6.6 x 10 20 m- 3 . Ver- 
tical solid line indicates the scaled enhancement frequency 
zuo — 1.22, where e" oc (w) vanishes, (a) e = ei oc (ci7o) = 5.2 
and (b) e = £i oc (ci7o) = 8.0. 

vortices distributed periodically and an envelope density 
profile which decreases as moving away from the center 
of the trap. The envelope is a slowly varying function 
compared to the periodicity of the vortices. Radius of 
the cloud is much greater than the periodicity, such that 
there may be few hundred vortices which are experimen- 
tally observable. Moreover, the distortion of the lattice 
near the edges of the condensate is small. 

In our past work, Q , we have numerically investigated 
the effects of the finite size and imperfections in the pe- 
riodicity. We observed that the positions of the gaps 
are not strongly effected, despite the occurrence of extra 
scattering due to the smooth density envelope over the 
lattice. 

Thus, in this paper we consider an infinite homoge- 
nous vortex lattice and concentrate on the effects of the 
frequency dependence of the dielectric function. 

A two-dimensional photonic crystal supports only two 
polarization modes for the in plane propagation of light 
[l~2l |. If the magnetic field, H, is perpendicular to the 
plane of periodicity, this mode is called transverse electric 
(TE). In a TE mode electric field, E, is perpendicular to 
the axis of vortices. Similarly, the mode with E parallel 
to the vortex axis is called transverse magnetic mode 
(TM). 

Let's first focus on the TE modes. We take the vortices 
to be aligned in the z direction, forming a periodic array 
in the x-y plane. A generalized eigenvalue equation for 

H || z 

Vx(^Vx^)=0^(r1 (8) 

is derived by decoupling the Maxwell equations for H, af- 
ter the substitution D(r,u) = e(r,u>)E(r,ui) [l5j|. Unlike 
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the frequency independent case, differential operators on 
the left hand side also depend on uj. Moreover, differen- 
tial operator is not Hermitian because of the imaginary 
part of the dielectric function ([5]). This causes eigenfre- 
quencies to be complex. Since, in general, Eq. ([8]) is not 
analytically solvable, we determine the eigenfrequencies 
computationally by plane wave expansion. 

Using Bloch-Floquet Theorem [25| the magnetic field 
can be expressed in terms of the reciprocal lattice vectors 
G as 



Similarly, inverse dielectric function is expanded as 
f 



e(r, uj) 



= J2 £ G' e[G ' p 



(9) 



(10) 



where the Fourier components e^, are 



£ G'M = 1 



i r e 



-iG'-r 



A J e(r,w) 



dV 



(11) 



The integration is carried out over the Wigner-Seitz unit 
cell of area A. 

We substitute the expansions ([9]) and (|10[) into the 
Master equation j8]), and obtain the expression 



C (G- 

G 



*_ &) {w)H d , [(k + G)-(k + G')] = (^y 



H G 



(12) 

We note that, if the dielectric function were real, e% - = 

G — G 

£ g'-G' ^ ne ma trix, represented by equation p2[) . would 
be real. Eigenfrequencies would also be real. Moreover, 
due to inversion symmetry of the unit cell eg's would be 
real. However, the presence of complex dielectric func- 
tion destroys the Hermiticity. Eigenfrequencies are, in 
general, complex. 

Real space basis vectors for a triangular lattice are 

ai = ax and c?2 = ai^x—^y). Corresponding re- 
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ciprocal lattice basis vectors are b\ = fco ^ 

and b 2 = hoy, where magnitude of both vectors are 
fco = (2/\/3)(27r/a). Any lattice point, in the summa- 
tion, can be written as G — n\bi + 77,262, where n\ and 
ri2 are all integers. We also denote the fourier compo- 
nents of inverse dielectric function as eg = £ ni ,n 2 - 

For computational purposes, we limit the number of 
G vectors over which the summation will be carried out. 
We consider a parallelogram in the reciprocal space over 
which ni and n 2 runs through —N to N. TV is a positive 
integer. This gives a (27V + l) 2 x (2N + l) 2 matrix of 
elements 

Mij{oj) = s v1:V2 (cj) (k + nifei + n 2 b 2 ^ ■ (k + n'-J)i + n' 2 6 2 ) 



where the dependence of indices are given by 

i = (2N +\)n 1 +n 2 and j = (2N + l)n[ + n' 2 . (14) 
We use the notations 



11 — rii and r/2 



n 2 



(15) 



Solution of Master equation ([5]) reduces to the deter- 
mination of eigenfrequencies u> for each wave vector k. 
All distinct values of uj are obtained by choosing fc in the 
first Brillouin zone. 

When the dielectric function is independent of to, the 
eigenfrequencies would be easily determined by straight- 
forward matrix diagonalization [l2| . However, the depen- 
dence of Fourier elements eq, (uj) on the frequency forces 
the calculations to be carried out by relying on the con- 
dition of vanishing determinant, 



det(M) = 



(16) 



(13) 



Numerical calculations are based on finding the zeros of 
the determinant of matrix M, as a function of fc or lu. 
The zeros of the complex function is computed using a 
least squares method. We have checked the convergence 
of the solutions using different initial points. 

In the constant dielectric case one chooses a fc value as 
the input and determines the u> value. This is because 
uj is only on diagonals while fc is in every element of the 
matrix. However, the situation is completely different 
in the frequency dependent case. Both fc and uj exist in 
every element of the matrix. One may solve fc for the 
input values of uj, as well as determining the uj values 
entering the fc as input. Two methods reveal different 
physical pictures [15| . 

Choosing real u> values as input, one determines, in 
general, complex fc values whose imaginary part gives 
the spatial attenuation of the propagating wave. On the 
other hand, entering real fc values one solves for complex 
u, whose imaginary part determines the temporal atten- 
uation of the wave. Since we are mainly interested in 
the spatial attenuation of the waves, we followed the first 
method. However, we checked that the two approaches 
give parallel results. As a result of this procedure we ob- 
tain complex wave vector values. We denote the real and 
imaginary parts of the wave vector as fc = ku + iki- 

We note that fc/ value may imply two different phe- 
nomena: reflection or absorption. If the dielectric func- 
tion is real, imaginary part of the wavevector, fc/, has 
a simple interpretation. The incident wave is totally re- 
flected while penetrating into the crystal up to a distance 
of 27r/fc/. However, if dielectric function is complex one 
cannot distinguish between reflection and absorption for 
a given value of fc/. A mixture of both occurs. 

In our computations, we use 11 x 11, parallelogram 
shaped grid of plane waves. This corresponds to a 121 x 
121 dimensional matrix of elements given in Eq. (|13p . 
We define the determinant of the matrix as a function of 
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k and solve for the zeros of this complex function. We 
determined the complex k values corresponding to real 
u>. Resulting band structures are plotted in Figs. [4] and 

El 

Although we described our method for the TE modes, 
TM modes can be calculated similarly. 




IV. RESULTS AND DISCUSSION 

When the dielectric function is frequency independent, 
the Master equation © is scalable. That is, the struc- 
ture of photonic bands, expressed in terms of scaled fre- 
quency u>' = uja/2irc, is independent of dimensions of the 
unit cell. It only depends on the lattice type and the 
internal structure within the unit cell. In the case of a 
frequency dependent dielectric, however, such a scaling is 
not possible, as a new length scale is introduced into the 
system. Thus, we first calculate the band structure with 
constant e, and then discuss the change due to frequency 
dependence of the dielectric constant. 

We give the constant dielectric photonic bands of a 
triangular lattice of Rubidium gas for two different sets 
of parameters, Fig. [3] The first set, N = 5.5 x 10 20 m~ 3 
and a — 10£, has a directional pseudo-gap in the TM 
direction, Fig. EK, when the dielectric constant is chosen 
as its value at the enhancement frequency, e = e loc (-nc7o) = 
5.2. The pseudo-gap lies between the frequencies u> = 
0.27-0.31 (27rc/a), with its center at uj 9 = 0.285(27rc/a). 

The second set of parameters, N — 6.6 x 10 20 m -3 
and a = 4.5£, is chosen so that there is a complete pho- 
tonic band gap when the dielectric constant is at its en- 
hanced value e = e loc (n7o) = 8.0. The band gap lies 
between u> — 0.30-0.32(27rc/a) with mid gap frequency 
LJg = 0.31(27rc/o). 

Strong frequency dependence of the dielectric suscep- 
tibility, Fig. El will modify the structure of the bands 
significantly. We note that, the dielectric function © 
is different from one only in the frequency range of 
u> = Qq ± 57. The natural lattice frequency 2irc/a that 
we used in the scaling of Fig. El is 7 orders of magni- 
tude greater than the decay rate, 7. For typical values 
of the lattice parameter in a rotating BEC, a ~ 200 nm, 
lattice frequency comes is 2nc/a ~ 10 15 Hz, whereas the 
decay rate is only 7 = 2n x 6 x 10 6 Hz. The bands 
of frequency dependent e loo (u>) will be different than the 
propagation in vacuum, for only about ~ IO7 around the 
enhancement frequency fio- On the other hand, index 
enhancement without absorption is achievable in a more 
narrow range of frequency, about O.I7. 

Still, the constant e bands give us an idea about how 
to arrange the lattice parameter, a, to obtain a band gap 
with the frequency dependent e Iot .(u>). In order to obtain 
a gap, we must arrange the enhancement frequency fio 
such that, it lies in the band gap of the corresponding 
constant dielectric case. A good choice is to place f2o at 
the center of the band gap, u> g . Thus, we tune the lattice 




FIG. 3: TE modes of a triangular vortex lattice with fre- 
quency independent e. Dielectric constants and lattice pa- 
rameters are (a) e = 5.2 and a = 10£, (b) e = 8 and a = 4.5£. 
Filling fractions of vortices, / = (2tt/V3) x {R 2 /^ 2 ) with effec- 
tive radius R ~ 2£, are 15% and 71%, respectively. Dielectric 
constant is the value of dielectric function ([6j| at the enhance- 
ment frequency, e = e loc (-cuo). Density profile of the unit cell 
is treated using the Pade approximation [24|. (a) There ex- 
ists a a directional pseudo-band gap with midgap frequency 
at u' g = 0.285. (b) There is a complete band gap with gap 
center at uo' g — 0.31. 



parameter a such that f2o = u) g (2Ttc/a), which gives 

27TC 



a = uj 



9 LUab + ^07 



(17) 



were u' g is obtained from constant dielectric calculations 
as in Fig. El 

In conventional photonic crystals it is generally not 
possible to change the lattice parameter, once the sample 
is manufactured. However, in the case of a rotating BEC, 
spacing between the vortex cores is continuously tunable. 
Density of the vortices depends on the rotation frequency, 
so the lattice parameter a can be decreased/increased by 
increasing/decreasing the rotation rate. Filling factor / 
of the lattice depends on / ~ (£/a) 2 , as coherence length 
£ determines the vortex core radius (see Eq. |7])). Co- 
herence length £ can be adjusted by density N. Alter- 
natively, £ can be adjusted by controlling a sc via Fesh- 
bach resonances [26j]. We note that, one might be able 
to design more convenient sets of parameters for specific 
experiments. Using different alkali atoms, like cesium, 
stronger index contrasts can be achieved due to larger 
transition wavelengths. By employing different index en- 
hancement schemes, such as the Raman scheme, broader 
index enhancement windows could be translated to wider 
band gaps. 

Choosing lattice parameter a as in Eq. (fT?)) . we cal- 
culate the photonic bands for the frequency dependent 
dielectric function ([5]). Away from the enhancement 
frequency tvq, the dielectric function is complex, (jll2p . 
Band structures are depicted in figures |H and E] for the 



6 



same density and filling factor parameters with figures^ 
andEJa, respectively. In both figures |H and El the real and 
the imaginary parts of the wave vector, ku and ki, are 
displayed separately. Enhancement frequency vjq = 1.22 
is marked in all plots. The lattice parameters a are cho- 
sen as a — 226 nm in Fig. [4] and a = 246 nm in Fig. 

El 

For a real dielectric function, it is very easy to identify 
the band gaps. The wave vector k is real when there is 
propagation, and complex (with ku on the band edge) 
if the frequency is in a band gap. For a complex e(u>), 
however, identification of band gaps is not straightfor- 
ward. One can determine the existence of a band gap 
by considering the frequency values where e is real. If a 
nonzero ki is present, then there exists a band gap at that 
frequency. However, the width of the band gap cannot 
be directly identified by considering only the k\ values 
away from the enhancement frequency. For a complex e, 
nonzero value of ki may be due to the absorption as well 
as the effect of the band gap. Thus, we first discuss the 
existence of the band gaps in Fig. |4]& Fig. [5] and discuss 
the gap widths later. 

In Fig. [H), imaginary parts of the wave vector kj is 
plotted for different propagation directions. In the TM 
and TK directions only a single band exists within the 
enhancement window, while for the MK direction there 
are two bands. At the enhancement frequency vjq = 1.22, 
two of these bands have zero kj, while the other two have 
a complex wavevector. In accordance with the discussion 
in the previous paragraph, we identify the existence of a 
pseudo-band gap in the TM propagation direction. Thus, 
incident light (exactly at tjjq) would propagate in the 
TK, MK directions while it would be stopped in the TM 
direction. 

In the second case, Fig. El all of the four bands have 
nonzero kj at vj = vjq. This indicates the existence 
of a complete band gap at the enhancement frequency. 
Incident light is stopped for all propagation directions. 

We see that the conclusions for the existence of pho- 
tonic band gaps obtained by constant e calculations are 
not modified, even when the strong frequency depen- 
dence of e(cj) is taken into account. 

Although the existence of directional and complete 
band gaps are demonstrated at the enhancement fre- 
quency, the widths of these gaps cannot be determined 
by only investigating the behavior of k\. In the vicinity 
of vj = tuq , one cannot distinguish wether the wave is de- 
caying (e~ kl ' r ) due to absorption, or due to the existence 
of a band gap. 

To be able to define the width of the gap, we calculate 
the behavior of the Poynting vector 

S(r) = l -E{r) x H(f)* (18) 

in the crystal. Real part of the Poynting vector, Sn(r), 
gives the energy flux of the field at position r. Imaginary 
part Si(f) is a measure of the reactive (stored) energy 
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FIG. 4: (a) TE modes of triangular vortex lattice with fre- 
quency dependent dielectric function e loc (vj) (Fig. |2j, and (b) 
imaginary parts of the wave vector ki corresponding to each 
mode. Particle density is TV = 5.5 x 10 20 m -3 and lattice con- 
stant is a = 10£. Enhancement frequency f2o is tuned to the 
band gap at the M edge (uo g — 0.285(27rc/a)) of the constant 
dielectric case (Fig. [3^0 . MK bands are plotted in a limited 
region, because of high ki values out of the given frequency 
region. There exists a directional gap in the TM propagation 
direction. 

[27| . For a frequency dependent, but real, e(uj) Poynt- 
ing vector is pure imaginary in the band gaps and real 
otherwise. However, for complex e(uj) imaginary part of 
the S may also be due to absorption. Although we make 
the similar statements on ki and Si, together they are 
sufficient to determine the width of the band gap. 

Using the results of our band structure calculations, 
we compute the average of the Poynting vector, (S), in 
the unit cell. We define 

a = \{S z )\/\(S)\ , (19) 

which corresponds to the rate of reactive energy. Aver- 
ages of (Si) and (S) are computed along the TM direction 
in order to investigate the gap width in Fig. [5] 

Fig. [6]displays a marked increase in the reactive energy 
ratio near the enhancement frequency. This increase can 
not be directly caused by the imaginary part of e IOC (tzj), 
as near the enhancement frequency this imaginary part is 
decreasing to zero. Thus, the peak in the reactive energy 
ratio must be caused mainly by the periodicity of the 
crystal. Presence of the band gap increases the reactive 
energy ratio despite decreasing absorption. 

We define the width of the photonic band gap as the 
full width at half maximum of the peak at the enhance- 
ment frequency. With this definition we find an effective 
gap in the frequency range w — zuq ± 0.043. In familiar 
units this translates to a band width of 3.30 MHz. Using 
the same method we find a band gap of width 5.98 MHz 
for the parameters of Fig. El 
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FIG. 5: (a) TE bands of triangular vortex lattice with fre- 
quency dependent dielectric function ei oc (ro) (Fig. [2p), and 
(b) imaginary parts of the wave vector ki corresponding to 
each mode. Particle density is N = 6.6 x 10 20 m -3 and lattice 
constant is a — 4.5f. Enhancement frequency Qo is tuned to 
the band gap at the M edge (ui g — 0.31(27rc/a)) of constant 
dielectric case (Fig. [3p). There exists a complete band gap. 




We also performed similar calculations for the TM 
modes and obtained similar band structures. TM modes, 
also, give directional and complete band gaps when the 
lattice parameter is properly tuned. However, since the 
band gaps of TE and TM modes do not coincide in gen- 
eral, one can not obtain band gaps for both modes, with- 
out further tuning. 

V. CONCLUSION 

We calculated the photonic bands for an index en- 
hanced vortex lattice, considering a frequency dependent 
dielectric function. Our motivation was the possibility 
of the direct measurement of the rotation frequency in 
Bose-Einstein condensates using the directional band gap 
of the photonic crystal. We validated the main conclu- 
sion of our previous work that photonic band gaps 
can be created via index enhancement on vortex lat- 
tices of BECs. Specifically, we presented two examples 
showing that both directional and complete band gaps 
are possible within experimentally realizable parameter 
regimes. For the specific parameters and index enhance- 
ment scheme we considered, band gaps of order a few 
MHz width are obtained. We also discussed how band 
gaps are designed for specific parameter values, and how 
band gap widths can be increased. 

Unlike the previous calculations performed for metal- 
lic photonic crystals, complex dielectric function varies 
rapidly with the frequency for index enhanced media. 
Strong frequency dependence is due to the high dipole 
moment, established trough atomic coherence in a nar- 
row frequency range. To our knowledge, such a periodic 
structure composed of index-enhanced media, is investi- 
gated for the first time. We showed that the photonic 
band structure in such a medium can be reliably cal- 
culated by numerically computing the zeros of the de- 
terminant of the Master equation. We also developed a 
method, based on the calculation of Pointing vector, to 
determine the effective widths of the photonic band gaps 
in media with frequency dependent complex dielectric 
function. 



FIG. 6: Reactive energy ratio a for the TM band of Fig. 
[3] Vertical dashed line indicates the enhancement frequency 
zuo = 1.22. Shaded region is the effective photonic band gap. 
Width of the peak determines the width of the gap to be 
ill = fl ± 0.0437 which corresponds to ±1.65 MHz. 
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